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I review some theoretical aspects of neutrino oscillations in the case when more than two neutrino flavours are 
involved. These include: approximate analytic solutions for 3-fiavour (3f) oscillations in matter; matter effects 
in u,j, «-» v T oscillations; 3f effects in oscillations of solar, atmospheric, reactor and supernova neutrinos and in 
accelerator long-baseline experiments; CP and T violation in neutrino oscillations in vacuum and in matter; the 
problem of C/ e 3; 4f oscillations. 



1. INTRODUCTION 

Explanation of the solar and atmospheric neu- 
trino data in terms of neutrino oscillations re- 
quires at least three neutrino species, and in fact 
three neutrino species are known to exist - v e , 
and v T , If the LSND experiment is correct, then 
probably a fourth neutrino type should exist, a 
light sterile neutrino v s . However, until relatively 
recently most of the studies of neutrino oscilla- 
tions were performed in the 2-flavour framework. 
There were essentially two reasons for that: (1) 
simplicity - there are much fewer parameters in 
the 2-flavour case than in the 3-flavour one, and 
the expressions for the transition probabilities are 
much simpler and by far more tractable, and (2) 
the hierarchy of Am 2 values, which allows to ef- 
fectively decouple different oscillation channels. 
The 2-flavour approach proved to be a good first 
approximation, which is a consequence of the hi- 
erarchy Atuq <C Am atm and of the smallness of 
the leptonic mixing parameter \U e s\. 

However, the increased accuracy of the avail- 
able and especially forthcoming neutrino data 
makes it very important to take into account even 
relatively small effects in neutrino oscillations. In 
addition, the experimentally favoured solution of 
the solar neutrino problem is at present the LMA 



MSW one, which requires the hierarchy between 
Atoq and Am atm to be relatively mild. Also, ef- 
fects specific to > 3 flavour neutrino oscillations, 
such as CP and T violation, are now being very 
widely discussed. All this makes 3-flavour (or 4- 
flavour) analyses of neutrino oscillations manda- 
tory. 

In my talk I review some theoretical issues per- 
taining to neutrino oscillations in the case when 
more than two neutrino species are involved. I 
mainly concentrate on 3-flavour (3f) oscillations 
and only very briefly consider the 4f case. The 
topics that are discussed include: approximate 
analytic solutions for 3f oscillations in matter; 
matter effects in za, <-> v T oscillations; 3f ef- 
fects in oscillations of solar, atmospheric, reactor 
and supernova neutrinos and in accelerator long- 
baseline experiments; CP and T violation in neu- 
trino oscillations in vacuum and in matter; the 
problem of U e 3\ 4f oscillations. 

2. NEUTRINO OSCILLATIONS IN 
MATTER (3f) 

Neutrino oscillations in matter are described 
by the evolution equation i(d/dt)v — Hu, where 
v = (y e v T ) T and 



*On leave from National Research Centre Kurchatov In- 
stitute, Moscow, Russia. Supported by the Calouste Gul- 
benkian Foundation as a Gulbenkian Visiting Professor. 



H = 



E x 

r\ o e 2 o 

E 3 



V{t) 





(1) 



2 




Figure 1. Normal mass hierarchy 




Figure 2. Inverted mass hierarchy 



The effective potential V — \/2GpN e is due to 
the charged-current interaction of v e with the 
electrons of the medium. The neutral current in- 
duced potentials are omitted from Eq. ([!]) be- 
cause they are the same for neutrinos of all three 
species and therefore do not affect neutrino os- 
cillations. This, however, is only true in leading 
(tree) order; radiative corrections induce tiny dif- 
ferences between the neutral current potentials of 
v e , and v r and, in particular, result in a very 
small - v T potential difference V^r ~ 10~ 5 V 
|l| . This quantity is negligible in most situations 
but may be important for supernova neutrinos. 

For matter of constant density, closed-form so- 
lutions of the evolution equation can be found 
however, the corresponding expressions are 
rather complicated and not easily tractable. For 



a general electron density profile N e ^ const no 
closed-form solutions exist. It is therefore desir- 
able to have approximate analytic solutions of the 
neutrino evolution equation. A number of such 
solutions were found, most of them based on the 
expansions in one (or both) of the two small pa- 
rameters: 

Arr4/Arr4 = Am%/Aml tm < 0.1 , (2) 

|[/ e3 | = |sin0 13 | <0.2 [3]. (3) 

Our numbering of neutrino mass eigenstates cor- 
responds to that in Figs. [l] and g, which also 
show schematically the possible neutrino mass hi- 
erarchies and the flavour composition of neutrino 
mass eigenstates. 

In the limits Am^ = or U13 = the tran- 
sition probabilities acquire an effective 2f form. 
When both these parameters vanish, the genuine 
2f case is recovered. 

2.1. Constant-density matter 

In the case of matter of constant density ap- 
proximate solutions of the neutrino evolution 
equation were found using the expansion in a = 
A/rig/Am^tm in An expansion in both a and 
sin #13 was used in pj. The 24 <-> transition 
probability found in [g| has the general form 

P(v e «-» v^) ~ sl 3 P 2 (Am 2 31 ,e 13 ,N e ) + 

c^fAm^^.JVe) + interf. term, (4) 

where the quantities Pi are the 2f transition prob- 
abilities in matter depending on the correspond- 
ing parameters shown in the parentheses. The 
interference term, which is linear in both a and 
sin6>i3, describes the genuine 3f effects, both CP- 
conserving and CP-violating. 

2.2. Arbitrary density profile 

Matter of constant density is a good first ap- 
proximation for long-baseline accelerator neu- 
trino experiments (neutrinos traverse the man- 
tle of the Earth). However, it is not very use- 
ful for describing the oscillations of solar, atmo- 
spheric and supernova neutrinos. An alternative 
approach is to consider matter with an arbitrary 
density profile and reduce the problem to an ef- 
fective 2f one plus easily calculable 3f corrections. 
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Figure 3. LMA allowed parameter region for 

013 = @ 

This has been done using the expansion in a in || 
and the expansion in sin #13 in |Q-fj. A different 
approach, based on the adiabatic approximation, 
was employed, e.g., in fio|| . 

2.3. Matter effects in v T oscillations 

Since the matter-induced potentials for and 
v r are the same (neglecting the radiative correc- 
tions), in the 2f case the «-> v T oscillations are 
not affected by matter. This, however, is not true 
in the 3f case; therefore matter effects on «-> v T 
oscillations is a pure 3f effect. It vanishes only 
when both Am^i and U e3 vanish. 

3. 3f EFFECTS IN v OSCILLATIONS 

We shall now discuss 3f effects in oscillations of 
neutrinos from various sources. 

3.1. Solar neutrinos 

In the 3f case, solar v e can in principle oscillate 
into either v^, or v T , or some their combination. 
What do they actually oscillate to? 

It is easy to answer this question. The small- 
ness of the mixing parameter \U e 3\ implies that 
the mass eigenstate v 3 is approximately given by 
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Figure 4. LMA allowed parameter region for 
sin 2 6»i3 = 0.04 

IAj ~ S 2 3 + c 23 v T (5) 

and, to first approximation, does not participate 
in the solar neutrino oscillations. From the uni- 
tarity of the leptonic mixing matrix it then follows 
that the solar neutrino oscillations are the oscilla- 
tions between v e and a state v' which is the linear 
combination of and v r , orthogonal to v 3 : 

v' = c 23 ffi - S23 v T (6) 

Since the mixing angle 6*23, responsible for the 
atmospheric neutrino oscillations, is known to be 
close to 45°, Eq. (g) implies that the solar v e 
oscillate into a superposition of and v T with 
equal or almost equal weights. 

What are the 3f effects in the oscillation prob- 
abilities? Since at low energies and v T are 
experimentally indistinguishable, all the observ- 
ables depend on just one probability - the v e sur- 
vival probability P{v e — ► v e ). Averaging over fast 
oscillations due to the large mass squared differ- 
ence Am 2 tm = Am§ x yields Q 

P{v e -» u e ) ~ ci 3 P 2ee {Am 2 21 ,6 12 ,N eg ) + sj 3 . (7) 

Here P2ee(Am2 1 , 9i2, N c s) is the 2f survival prob- 
ability of v e in matter with the effective electron 
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density N cS = cf 3 N e . 

As follows from the CHOOZ data f|, the sec- 
ond term in Eq. (Q), sf 3 , does not exceed 10~ 3 , 
i.e. is negligible. At the same time, the coeffi- 
cient cf 3 of P 2ee in the first term may differ from 
unity by as much as ~ 5 - 10 %. Thus, 3f effects 
may lead to an energy-independent suppression of 
the v e survival probability by up to 10%. With 
high precision solar data this must be taken into 
account. This is illustrated by Figs. || and || 
fL2| ): The difference between the cases 6\ 3 = 
and sin 2 # 13 = 0.04 (which is about the maximum 
allowed by CHOOZ value) is quite noticeable. 

3.2. Atmospheric neutrinos 

(1) The dominant channel <-> v T . In the 2f 
limit, there are no matter effects in this channel 
(neglecting tiny V^ r caused by radiative correc- 
tions). The oscillation probability is independent 
from the sign of Am^, i.e. cannot differentiate 
between the normal and inverted neutrino mass 
hierarchies. The 3f effects result in a weak sensi- 
tivity to matter effects and to the sign of Am 31 . 

(2) The subdominant channels v e <-» ^V,t- Con- 
tributions of these oscillation channels to the 
number of fi - like events are subleading and 
difficult to observe. For e-like events, one could a 
priori expect significant oscillations effects. How- 
ever, these effects are in fact strongly suppressed 
because of the specific composition of the atmo- 
spheric neutrino flux and proximity of the mixing 
angle 6*23 to 45°. Indeed, in the 2f limits one finds 



F2 



P 2 (Am 2 31 ,e 13 ,V)-(r S 2 23 -l) (8) 



in the limit Ar 



and 



F - F Q 
F9 



= P 2 (Am 2 21 ,9 12 ,V)-(rc 2 23 -l) (9) 



in the limit s± 3 — > [Q. Here F® and F e are 
the v e fluxes in the absence and in the presence 
of the oscillations, respectively, and r = F®/F® . 
At low energies r ~ 2; also, we know that s 23 ~ 
c 23 ~ 1/2. Therefore the factors (rs 23 — 1) and 
( rc 23 — 1) in Eqs. (||) and (||) are very small and 
strongly suppress the oscillation effects even if the 
transition probabilities P 2 are close to unity. This 



happens because of the strong cancellations of the 
transitions from and to the v e state. 

All this looks as a conspiracy to hide the oscil- 
lation effects on the e-like events! This conspir- 
acy is, however, broken by the 3f effects. Keeping 
both Aniji and si 3 in leading order yields || 

^-fl ~ P 2 (Aml,e 13 )-{rs 2 23 -l) (10) 

+ P 2 (Am 2 21 ,d 12 )-(rc 2 23 -l) (11) 
- 2s 13 s 23 c 23 rRe(A* e A fie ) 

The interference term, which represents the gen- 
uinely 3f effects, is not suppressed by the flavour 
composition of the atmospheric neutrino flux; it 
may be responsible (at least, partially) for some 
excess of the upward-going sub-GeV e-like events 
observed at Super-Kamiokande |J. 

3.3. Reactor antineutrinos 

Since the average energy of reactor i> e 's is E ~ 3 
MeV, for intermediate-baseline experiments, such 
as CHOOZ and Palo Verde (L ~ 1 km), one has 
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This justifies the use of the one mass scale domi- 
nance approximation, which gives 



P{v e — * u e ) = 1 - sin 2 2^i3 • sin 2 
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L ,(13) 



a pure 2f result. However, in the case of the LMA 
solution of the solar neutrino problem, at high 
enough confidence level Am 21 can be comparable 
with Am 31 , and the second condition in (12) may 
not be valid. In such a situation the 3f effects 
coming through the subdominant Am 21 should 
be taken into account, The analyses jl3| show 
that the constraints on |J7 e 3| derived from the 
CHOOZ experiment become slightly more strin- 
gent in that case. However, the new SNO data 
PI disfavour large values of Am 21 and so make 
this possibility less likely. 

For KamLAND, which is a very long baseline 
reactor experiment (L ~ 170 km), one has 
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L > 1 (for LMA) . (14) 
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Averaging over the fast oscillations driven by 

A "4m y ields 



P(y e -» v e ) = c? 3 P 2gg (Am^,0 12 ) + s? 3 . (15) 

This has the same form as Eq. (|7|), except that 
the 2f survival probability P2ge has to be calcu- 
lated in vacuum rather than in matter; it is in fact 
given by Eq. ( |l3| ) . The probability ( |l5| ) can differ 
from the 2f probability (ph by up to - 10%. 



3.4. LBL accelerator experiments 

(1) Up disappearance. 

3f effects can result in up to ~ 10% corrections 
to the disappearance probability mainly due to 
the factor cf 3 in the effective amplitude of the 



Vt oscillations, 



sin 2 (26» MT ) 



13 



sin 2 29 
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(16) 



Another manifestation of 3-flavourness are small 
matter effects in <-> v T oscillations. The same 
applies to v T appearance in experiments with the 
conventional neutrino beams, disappearance 
also receives contributions from the subdominant 
Up <-> v e oscillations. 

(2) Up appearance at neutrino factories; u e ap- 
pearance at neutrino factories and in experiments 
with the conventional neutrino beams. 

These are driven by the v & «-> !/„ T oscillations. 
There are two channels through which these 
subdominant oscillations can proceed - those 
governed by the parameters (#13, Am^) and 
(012, Aniji). For typical energies of the LBL ac- 
celerator experiments (a few GeV to tens of GeV), 
and assuming the LMA solution of the solar neu- 
trino problem, one finds that for 6*13 in the range 
3 • 10~ 3 < 6<i3 < 3 • 10~ 2 the two channels com- 
pete; otherwise one of them dominates. 

Unlike in the case of atmospheric neutrinos, 
there is no suppression of the oscillation effects 
on the v e flux due to the flavour composition of 
the original flux. 

The dependence of the oscillation probabilities 
on the CP-violating phase <5cp (both ~ sin<5cp 
and ~ cos<$cp) comes from the interference terms 
and is a pure 3f effect. The 3f effects will be 
especially important for the future experiments at 
neutrino factories which are designed for precision 
measurements of neutrino parameters. 



3.5. Supernova neutrinos 

In supernovae, matter density varies in a very 
wide range, and the conditions for the three MSW 
resonances are satisfied (taking into account that 
due to radiative corrections V^ r 7^ 0), see Fig. ||. 
The hierarchy Am|x <C Am 31 leads to the ap- 
proximate factorization of transition dynamics at 
the resonances, so that the transitions, to first 
approximation, are effectively 2f ones. However, 
the observable effects of the supernova neutrino 
oscillations depend on the transitions between all 
three neutrino species. 

The Earth matter effects on supernova neutri- 
nos can be used to measure |f/ e 3| to a very high 
accuracy (~ 10~ 3 ) and to determine the sign of 
Am 2 ! . 

The transitions due to the — v r potential 
difference V^ T caused by radiative corrections 
may have observable consequences if the origi- 
nally produced and v r fluxes are not exactly 
the same fll6l. 




Figure 5. Energy level crossing scheme for super- 
nova neutrinos 



4. CP AND T VIOLATION IN 
v OSCILLATIONS IN VACUUM 

The probability of v a — ► v\, oscillations in vac- 
uum is given by 



P(v a ,t -+Vb',t) 



s^ Ubie -i Ei (t-t o)u * ai 



(17) 
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In the general case of n flavours the leptonic mix- 
ing matrix U a i depends on (n— l)(n — 2)/2 Dirac- 
type CP-violating phases {<5cp}- If neutrinos are 
Majorana particles, there are n — 1 additional, so- 
called Majorana-type CP-violating phases. How- 
ever, they do not affect neutrino oscillations and 
therefore I shall not discuss them. 

Under CP transformation, neutrinos are re- 
placed by their antiparticles (v a ,b *V&), which 
is equivalent to the complex conjugation of U a f. 



CP : v a . b <- 



V a ,b 



Ki ({<*cp} - -{*cp}) • (18) 



Time reversal transformation interchanges the 
initial and final evolution times to and t in Eq. 
jl7|), i.e. corresponds to evolution "backwards 
in time" . As follows from Eq. ( [17| ) , the inter- 
change to ^ t is equivalent to the complex conju- 
gation of the exponential factors in the oscillation 
amplitude. Since the transition probability only 
depends on the modulus of the amplitude, this 
is equivalent to the complex conjugation of the 
factors Ubi and U*i, which in turn amounts to 
interchanging a ^ b. Thus, instead of evolution 
"backwards in time" one can consider evolution 
forward in time, but between the interchanged 
initial and final flavours: 

T : t ^ t <^ v a <-> v h 

Uai^Kt ({S CP } -> -{<5 CP }) . (19) 

Under the combined action of CP and T: 

CPT : v a , h <-> v aih h to ^ t (v a <-> v b ) 

=> P(y a -» u b ) -» P{9 b -» D a ) . (20) 

From CPT invariance it follows that CP violation 
implies T violation and vice versa. 

CP and T violation can be characterized by the 
probability differences 
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P(Va 
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ab — 



P(y a -> - P(y h v a ) . 



(21) 
(22) 



From CPT invariance it follows that the CP- and 
T- violating probability differences coincide, and 
that the survival probabilities have no CP asym- 
metry: 



CP and T violations are absent in the 2f case, so 
any observable violation of these symmetries in 
neutrino oscillations in vacuum would be a pure 
> 3f effect. 

In the 3f case, there is only one CP-violating 
Dirac-type phase 8cp and so only one CP-odd 
(and T-odd) probability difference: 

AP e CP = Apcp = Ap cp = APj (24) 



AP 



sin 



4si2 C12 S13 Ci 3 s 2 3 c 23 sin<5 C p x 
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It vanishes 

• when at least one Am 





or 90° 



A pCP _ A pT . 
^^ab — ^^ab ) 



AF CP = 



(23) 



• when at least one 9 

• when (5 C p = or 180° 

• in the averaging regime 

• in the limit L — ► (as L 3 ) 

Clearly, this quantity is very difficult to observe. 

5. CP AND T VIOLATIONS IN 
v OSCILLATIONS IN MATTER 

For neutrino oscillations in matter, CP trans- 
formation (substitution v a <-> j? a ~j implies not only 
complex conjugating the leptonic mixing matrix, 
but also flipping the sign of the matter-induced 
neutrino potentials: 

CP : U„i - U* ai ({<5 CP } -> -{5cp}) , 

V(r) -> -V(r). (25) 

It can be shown Q that in matter with an ar- 
bitrary density profile, as well as in vacuum, the 
action of time reversal on neutrino oscillations is 
equivalent to interchanging the initial and final 
neutrino flavours. It is also equivalent to complex 
conjugating U a i and replacing the matter density 
profile by the reverse one: 

T : U ai -» U* ai ({8 CV } -> -{^ C p}) , 

V{r) -> V{r) . (26) 

Here 

V{r) = V2G F N{r) , (27) 

N(r) being the reverse profile, i.e. the profile that 
corresponds to the interchanged positions of the 
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neutrino source and detector. In the case of sym- 
metric matter density profiles (e.g., matter of con- 
stant density), N(r) =N(r). 

An important point is that the very presence 
of matter (with unequal numbers of particles and 
antiparticles) violates C, CP and CPT, leading to 
CP violation in neutrino oscillations even in the 
absence of the fundamental CP-violating phases 
{<5cp}- This fake (extrinsic) CP violation may 
complicate the study of the fundamental (intrin- 
sic) one. 

5.1. CP violation in matter 

Unlike in vacuum, CP violation in neutrino os- 
cillations in matter exists even in the 2f case (in 
the case of three or more flavours, even when all 
{<5cp} = 0): 

P{Va -» Vb) ¥= P("a -> Vb) ■ (28) 

This is actually a well known fact - for example, 
the MSW effect can enhance the v e <-» v„ oscilla- 
tions and suppress the 9 e <-» ones or vice versa. 
Moreover, in matter the survival probabilities are 
not CP-invariant: 



P(v a -> v a ) ^ P(v a -> v a ) . 



(29) 



To disentangle fundamental CP violation from 
the matter induced one in the LBL experiments 
one would need to measure the energy depen- 
dence of the oscillated signal or the signals at two 
baselines, which is a difficult task. The (difficult) 
alternatives are: 

• LBL experiments at relatively low energies 
and moderate baselines (E ~ 0.1 - 1 GeV, L ~ 
100 - 1000 km) @. 

• Indirect measurements through 

(A) CP-even terms ~ cos <5 C p JlJ| ; 

(B) Area of leptonic unitarity triangle Jl9| . 
CP violation cannot be studied in the supernova 
neutrino experiments because of the experimental 
indistinguishability of low-energy and v r . 

5.2. T violation in matter 

Since CPT is not conserved in matter, CP and 
T violations are no longer directly connected (al- 
though some relations between them still exist 
|||2(J). Therefore T violation in neutrino oscilla- 
tion in matter deserves an independent study. Its 
characteristic features are: 



• Matter does not necessarily induce T viola- 
tion (only asymmetric matter with N(r) ^ N(r) 
does) . 

• There is no T violation (either fundamental or 
matter induced) in the 2f case. This is a simple 
consequence of unitarity. For example, for the 
(z/ e , j/ju) system one has 



Pe, 
P, 



Pue = 1 



(30) 



from which P ejl = P^. 

• In the 3f case there is only one T-odd proba- 
bility difference for z/'s (and one for P's), irrespec- 
tive of the matter density profile: 

AP? = APj T = AP r T e . (31) 

This is a consequence of 3f unitarity |^l| . 

The matter-induced T violation is an interest- 
ing, pure >3f matter effect, absent in symmetric 
matter (in particular, in constant-density mat- 
ter). It does not vanish in the regime of complete 
averaging of neutrino oscillations || . It may fake 
the fundamental T violation and complicate its 
study, i.e. the extraction of Sqp from the exper- 
iment. The matter- induced T violation vanishes 
when either U e3 = or Am^ = (i.e., in the 2f 
limits) and so is doubly suppressed by both these 
small parameters. This implies that the perturba- 
tion theory can be used to obtain analytic expres- 
sions for the T-odd probability differences. The 
general structure of these differences is 

(32) 



AP^ = sin S C p - Y + cos S C p ■ X . 



Here the first term (oc sin^cp) is due to the fun- 
damental T violation, whereas the second term 
is due to the matter-induced one. In the adi- 
abatic approximation one finds || X = J c g x 
(oscillating terms) , where 

sin(2(9i - 20 2 ) 



JeS = Sl2 C12 S13 C 13 S 2 3 C23 



(33) 



sin 20i2 

Here 0i and 02 are the mixing angles in matter in 
the (1-2) sector at the initial and final points of 
neutrino evolution, respectively; 0i — 02 is there- 
fore a measure of the asymmetry of the density 
profile. J ff has to be compared with the vacuum 
Jarlskog invariant 



J = S12 C12 Si 3 cj 3 S 2 3 C23 SU1 <5 C P ■ 



(34) 
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We see that the factor sin(26*i — 2#2)/sin2#i2 in 
Jeff plays the same role as the factor sin Sep in J. 

In an asymmetric matter, both fundamental 
and matter-induced T violations contribute to the 
T-odd probability differences AP^. This may 
hinder the experimental determination of the fun- 
damental CP- and T- violating phase <5cp . In par- 
ticular, in the accelerator LBL experiments one 
has to take into account that the Earth's den- 
sity profile is not perfectly spherically symmetric. 
To extract the fundamental T violation, strictly 
speaking one would need to measure 



v a * v b) Prcviyb * v a) ■ 



(35) 



where Pdir and P rcv correspond to the direct and 
reverse matter density profiles. (An interesting 
point is that even the survival probabilities P MAl 
and P TT can be used for that p|| . The v e survival 
probability P ee is an exception because in the 3f 
case it does not depend on <5cp This, 
however, is not true if v s is present [pi). 

In practical terms, it would certainly be diffi- 
cult to measure the quantity in (35): It would 
not be easy, for example, to move CERN to 
Gran Sasso and the Gran Sasso Laboratory to 
CERN. Fortunately, this is not actually neces- 
sary - matter-induced T violation due to imper- 
fect sphericity of the Earth's density distribution 
is very small. It cannot spoil the determination 
of 6 C p if the error in <5 C p is > 1% at 99% C.L. |]. 

Can we study T violation in neutrino oscilla- 
tions experimentally? Because of problems with 
the detection of e ± this seems to be difficult, 
but probably not impossible. To study matter- 
induced T violation would be a harder task. 
T-odd matter effects are expected to be negligi- 
ble in terrestrial experiments. They cannot be 
observed in the supernova neutrino oscillations 
because of the experimental indistinguishability 
of low - energy and v r . It could, however, af- 
fect the signal from ^GeV neutrinos produced in 
the annihilations of WIMPs inside the Sun 125(1 • 



6. A HYMN TO U, 
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The leptonic mixing parameter U e ^ plays a very 
special role in neutrino physics. It is of particular 
interest for a number of reasons. 



First, it is the least known of leptonic mixing 
parameters: while we have (relatively small) al- 
lowed ranges for the other two mixing param- 
eters, we only know an upper bound on \U&\. 
Its smallness, which looks strange in the light of 
the fact that the other two mixing parameters, 
012 and 6*23, are apparently large, remains essen- 
tially unexplained. (There are, however, some 
ideas which relate the smallness of \U&\ to that 



ofAm|/AmL @27|). 

The smallness of U e z is likely to be the bot- 
tleneck for studying the fundamental CP and T 
violation effects and matter-induced T violation 
in neutrino oscillations. The same applies to the 
determination of the sign of Am 2 ^ in future LBL 
experiments, which would allow us to discrimi- 
nate between the normal and inverted neutrino 
mass hierarchies. Therefore it would be vitally 
important to know how small U e 3 actually is. 

The parameter [7 e3 can be efficiently used to 
discriminate between various neutrino mass mod- 
els [^8[^9). It is one of the main parameters 
that drives the subdominant oscillations of atmo- 
spheric neutrinos and is important for their study. 
It also governs the Earth matter effects on super- 
nova neutrino oscillations. 

And finally, U e s apparently provides us with 
the only opportunity to see the "canonical" MSW 
effect. While matter effects can be important 
even in the case of large vacuum mixing an- 
gles, the most spectacular phenomenon, strong 
enhancement of mixing by matter, can only oc- 
cur if the vacuum mixing angle is small. From 
what we know now, it seems that the only small 
leptonic mixing parameter is C/ e 3- 

All this makes measuring U e 3 one of the most 
important problems in neutrino physics. 



7. 4f OSCILLATIONS 

If the LSND experiment is correct, the oscil- 
lations interpretation of the solar, atmospheric 
and accelerator neutrino data would require three 
distinct values of Am 2 : Am^ <g; Am 2 tm <C 
A^lsnd- This would imply the existence of at 
least four light neutrino species, v e , v^, v r and 
v s . A possible alternative is a strong CPT viola- 
tion in the neutrino sector, leading to inequalities 
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of Am 2 in the neutrino and antineutrino sectors 
30 1 ; I will not discuss this possibility here. 

In general, the 4f neutrino oscillations are de- 
scribed by 6 mixing angles 6>y, 3 Dirac-type CP- 
violating phases and 3 values of Amfj, i.e. are 
quite complicated. Fortunately, there is a sim- 
plification: The data admit only 2 classes of 4f 
schemes, the so-called (3+1) and (2+2) schemes. 
In the (3+1) schemes, three neutrino mass eigen- 
states are close to each other while the fourth one 
is separated from them by a large mass gap. This 
mass eigenstate is predominantly v s with small 
admixtures of the active neutrinos: 



v s + 0(e) • (v e , v^, v T ), e < 1 



(36) 



whereas v\ , and are the usual linear combi- 
nations of v ei and v T plus small (~ e) admix- 
tures of v s . In this scheme the amplitude of the 



oscillations at LSND is 



sin 2 20 LSND =4|£/ e4 £/ M4 | 2 ~ e 4 



(37) 



Strong upper bounds on \U e i\ and \U ^ from 
v e and disappearance experiments make it 
rather difficult to fit the LSND data in the (3+1) 
schemes Q. 

In the (2+2) schemes, there are two pairs 
of mass eigenstates with relatively small mass 
squared differences, Am Q and Am 2 tm , between 
the states within the pairs and a large separation 
(Am 2 SND ) between the two pairs. The state 
is predominantly in the pair responsible for the 
atmospheric neutrino oscillations, whereas v e is 
mainly in the pair responsible for Vq oscillations: 

fatm osc. : v^^v', vq osc. : v e «-> v" , 

where 



v ~ v T + v s + 0(e) • v e 
v" ~ v T + c 5 v s + 0(e) • 

The amplitude of the 
LSND is 



sin 2 26» LS nd ~ e 2 



(38) 



v P oscillations at 



(39) 



It is only of the second order in e and, unlike 
in the (3+1) case, the LSND data can be easily 
fitted. 



However, the (2+2) schemes suffer from a dif- 
ferent problem. The fractions of v s involved in 
the oscillations of atmospheric and solar neutri- 
nos must sum to unity in these schemes 133] : 



K^K>| 2 + \(v s W')f 



c? + st = 1 . 



(40) 



This sum rule is in conflict with the atmospheric 
and solar neutrino data. Indeed, the Super- 
Kamiokande atmospheric neutrino data lead to 
the upper limits sin 2 £ < 0.20 at 90% C.L. and 
sin 2 £ < 0.26 at 99% C.L. Q. At the same time, 
the (pre-SNO neutral current) solar neutrino data 
data imply sin 2 £ > 0.7 (90% C.L.); sin 2 £ > 0.48 
(99% C.L.) for the LMA solution of the solar neu- 
trino problem E3]. The recently published SNO 
neutral current data will probably strengthen 
this limit. Therefore, the (2+2) scenarios are also 
strongly disfavoured by the data. 

In the 4f case, there may be interesting mat- 
ter effects on neutrino oscillations |35|. CP vi- 
olation is potentially much richer than in the 3f 
case: there are several CP-violating observables, 
and large CP-odd effects are possible (in general, 
there is no suppression due to small Am Q ). Also 
large T violation (both fundamental and matter- 
induced) can occur. 

8. CONCLUSIONS 

3f effects in solar, atmospheric, reactor and su- 
pernova neutrino oscillations and in LBL accel- 
erator neutrino experiments may be quite impor- 
tant. They can lead to up to ~ 10% corrections to 
the oscillation probabilities and also to specific ef- 
fects, absent in the 2f case. The manifestations of 
> 3 flavours in neutrino oscillations include fun- 
damental CP violation and T violation, matter- 
induced T violation, matter effects in v T 
oscillations, and specific CP- and T-conserving in- 
terference terms in oscillation probabilities. The 
leptonic mixing parameter U e 3 plays a very spe- 
cial role and its study is of great interest. 

In the 4f case, large CP violation and (both fun- 
damental and matter-induced) T violation effects 
are possible. However, 4f scenarios are strongly 
disfavoured by the data. 
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